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Abstract 

Two-stream instability requires, essentially, two things to operate: a relative flow be- 
tween two fluids and some type of interaction between them. In this letter we provide 
the first demonstration that this mechanism may be active in a cosmological context. 
Building on a recently developed formalism for cosmological models with two, inter- 
penetrating fluids with a relative flow between them, we show that two-stream instabil- 
ity may be triggered during the transition from one fluid domination to the other. We 
also demonstrate that the cosmological expansion eventually shuts down the instability 
by driving to zero the relative flow and the coupling between the two fluids. 
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Two-stream instability of two, interpenetrating plasmas is a well-established phe- 
nomenon (11,3]. It is expected also to play a role in two, interpenetrating fluids, such as 



superfluid helium 0], gravitational-wave instabilities in neutron stars y]], and glitches 
in pulsars (due to the interior superfluid components) IslQl. It is not surprising that the 
instability shows up in diverse settings, since the basic requirements for it to operate 
are fairly generic: there must be a relative flow and some type of coupling between 
the two fluids. A window of instability may be opened when a perturbation of the flu- 
ids (hereafter, "mode") appears to be "left-moving", say, with respect to one fluid, but 
"right-moving" with respect to the other. The mechanism may be well-known, yet it 
was considered in a relativistic context only recently. Samuelsson et aJ. 1 1311 demon- 
strated quite generally the existence of (local) two-stream instability for a system of 
two, general relativistic fluids (using only causality - mode-speed less than one, in ge- 
ometric units - and absolute stability - modes are real for the fluids at rest - to constrain 
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the fluid properties). In this letter we consider a new application by demonstrating that 
two-stream instability may be triggered in cosmological settings. 

Although we use a specific example in this letter to exhibit the possibility of the 
mechanism to operate, we expect it to apply in various situations such as the presence 
of many inflaton scalar fields: the scalar fields turn into Bose-Einstein condensates 
behaving, for example, like the superfluids in a glitching neutron star, thus leading to 
two flows that have no particular reason to be aligned. Similarly, around the transition 
between dark matter and dark energy domination, the instability would provide a means 
of discriminating between a cosmological constant and a dark energy fluid. 

Observations over the past few decades have provided a wealth of information that 



can be used to constrain cosmological models III 2H . One of the more stringent is the 
leading-order observed homogeneity and isotropy of the Universe. Is it possible to have 
a relative flow in cosmologically important settings that would not have been detected 
yet? An ideal stage is during a cosmological transition between one phase of domi- 
nation to another. The requirement of a transition stems from the fact that one should 
naturally arrive at the current state of the Universe. We have shown, in a companion 
paper |9fl, that when one fluid flux dominates over the other (i.e. before and after the 
transition), one recovers the usual Friedmann-Lemaitre-Robertson- Walker (FLRW) be- 
havior. However, during the transition the model goes through a Bianchi I phase. We 
are thus led to ask whether the two-stream instability can be triggered during the same 
epoch? As we will show, the answer is "yes". As a proof-of-principle we consider a 
fairly simplified picture, leaving actual cosmological consequences and constraints on 
more elaborate models (that deserve further examination) for future work. 

Our cosmological two-fluid model has its relative flow along one direction, which 
we take to be the z-axis, cf. |@]. Orthogonal to the flow we impose two, mutually 
orthogonal spacelike Killing vector fields: one along the A-axis and another along the y- 
axis. The Killing vector fields X 1 * and V are thus X* =(0,1, 0, 0) and V = (0, 0,1,0). 
These two symmetries, and the remaining freedom in the choice of coordinates, imply 
the metric and fluid variables are functions of z and time t. The metric can therefore be 
written 

As 1 = -At 2 + A 2 .dx 2 + A 2 ,Ay 2 + A 2 dz\ ( 1 ) 

where the (still arbitrary) dimensionless functions A, (i e {x,y, z}) can in principle 
depend on both t and z- 

To demonstrate two-stream instability, it is sufficient to consider fixed "moments" 
of cosmological time {e.g. assume that the time-scale for the oscillations is much less 
than that of the cosmological expansion). For our mode analysis, this means any time 
derivatives of the background metric and the corresponding fluid flux will be ignored. 
We will also ignore the metric perturbations and relax the z-dependence in the back- 
ground metric and fluid flow so that our modes propagate in a spacetime of the well- 
known Bianchi I type. General cosmological perturbations in such a spacetime can be 



found in [14] 



As described in [9Q, we use the multi-fluid formalism originally due to Carter ]8[] 
(see J2l for a review). The fluid variables are two conserved four-currents, to be denoted 
n? - min (with g MV u^u v n = -1) and s M - su 1 ^ (with also g^u^u* = -1). The fluid 
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equations of motion are obtained from a Lagrangian A(n, s), which we take to be of the 
form of a fluid with constituent mass m coupled to a fluid with zero constituent mass; 
i.e. 

A(n, s) = -mn a - T n8 n' 7i - k s / , (2) 



where m, a, <r n , <x s , r ns , k s and /? are constants. The "bare" sound speeds II 1311 are given 
by 

2 = din/* 2 = dlnT 

n ~ fllnn* s ~ dins' ( ' 

where n = -dA/dn and T = -dA/ds are the associated conjugate momenta^ of the 
two fluids. The cross-constituent coupling reads 

_ 5 In// 7\v 
olni ju« 

Performing perturbations to linear order, the fluid densities and the z-component of 
the unit four-velocities take the form 

n(t,z) = nit) + 6Ne k - x4 ' , 

s(t,z) = sit) + SSe^*" , (5) 

where = ik t , 0, 0, k.) is the constant wave-vector for the modes and 5N, 6U^, etc. are 
the constant wave-amplitudes. Within the setting of Eq. (O, we see that the short- 
wavelength approximation for the modes is expressible as k t , k- » Hi = Aj/A,- for all 
i € {x,y,z}- 

The results in Figure Q] which displays the time evolution of the two flows M n s = 
A z u z ns /u' ns as well as the couplings (C n s,C sn ), show that the relevant coefficients are 
driven to zero by cosmological expansion, meaning, as expected and anticipated, that 
any instability will only operate during a finite time. Note that long before and after the 
transition, one fluid dominates over the other, so spacetime is effectively FLRW; only 
during the finite Bianchi I phase do the two fluids have comparable contributions 

Given a Bianchi I background, we place the terms of Eq. (0 into the Einstein 
and two-fluid equations of [91], expand, and keep only terms linear in the perturbed 
quantities, to arrive at the dispersion relation 



[iu n cr z - l) 2 c 2 n - icr z - u n ) 2 ] [iu s cr z - l) 2 c 2 - (cr z - m s ) 2 ] 

-C ns C sn iu n <r z - l) 2 iu s cr z - l) 2 = (6) 

for the mode speed cr z = -A z k,/k z . This relation is of the exact same mathematical 



form as the dispersion relation obtained by Samuelsson et al. 111311 [cf. their Eq. (69)]. 
It is thus immediately clear that our simplified model possesses all the ingredients for 
two-stream instability. 



1 We use notation that is reminiscent of the matter and entropy flows where n is the total matter number 
density, .v the entropy density, fi the chemical potential, and T the temperature. 



3 



The presence of the instability is shown in Figures [2] and [3] These figures provide 
graphs of the real and imaginary components of o~ z , versus the "e-folding" factor N 
defined as 

N = In (A x AyA z f /3 . (7) 

Figure [2] clearly shows four modes - denoted <r-(,) with i — 1, ... ,4 - all less than one, 
as should be the case. These modes evolve asymptotically towards constants corre- 
sponding to the "sound" speeds expected for both fluids. The overall asymmetry of the 
modes is due to the background relative flow. In the lower left corner of the figure, it 
is interesting to note the presence of a so-called avoided crossing, at which two of the 
modes "exchange" identity. 

The instability window is best seen in Figure [3] where the imaginary parts of <r z 
are graphed. Consistent with this is the "mode-merger" of the real parts that occurs in 
Figure|2] The results show that the window opens and closes during the same epoch as 
the transition takes place. Of course, the window closes because of the behavior shown 
in Figure[T] e.g. the relative velocity is quenched by the cosmological expansion. 

It should be clear that the mechanism discussed here is robust and should apply 
equally well, for example, to a system of coupled condensates. As stressed at the 
beginning, any system with two or more fluids, and a (non-gravitational) coupling be- 
tween them, could be subject to two-stream instability. What is perhaps unique in 
the cosmological context is that the instability shuts down automatically, without any 
fine-tuning, so long as there is overall expansion. 

The cosmological expansion thus provides a means of both initiating and ending 
the instability. This is mostly due to the cosmological principle, which states that the 
Universe should be, at almost all times, well described by a FLRW model. This im- 
plies that only finite epochs can exhibit a different behavior, in the case at hand that of 
a non-isotropic Bianchi phase. There remains to address the question of the origin of 
the relative flow. One might argue, for instance, that the inflation era, almost by con- 
struction, drives any primordial anisotropy to zero. During many field inflation itself, 
instabilities of this sort could be spontaneously initiated, even for a limited amount of 
time, and lead to new constraints. 

In the companion paper j9i], we suggested several scenarios that lead to a non- 
vanishing relative flow, some being also related to the question of the cosmologically 
coherent magnetic fields which are thought to exist 150 . These may actually reverse 
the question: it is conceivable that any model aimed at producing primordial magnetic 
fields on sufficiently large scales will also induce coherent anisotropies on these scales 
(the metric we used here is thus meant to describe only a cosmologically small region of 
space). The mechanism producing these magnetic fields should thus be tested against 
the instability suggested here. This question is quite natural given that two-stream 
instability is well-established for plasmas, and there is no reason why it should matter 
that these are placed in a cosmological setting. 

Based on the results presented here, we suggest that cosmological two-stream in- 
stability should be taken into account in further studies of all the transitions that could 
lead to its occurrence. In particular, the supposedly latest such transition, that ended the 
matter-domination era to the ongoing accelerated phase, could lead to drastically differ- 
ent observational predictions if the latter was driven by a mere cosmological constant 
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or by a cosmic fluid |@]: this new fluid would have no particular reason to be aligned 
along the matter flow, and hence an instability could develop, producing a character- 
istic anisotropy whose features still have to be investigated, at the typical scale corre- 
sponding to the transition. Some have argued that such an anisotropy has already been 
measured or that it could be using forthcoming Planck data lIToTl . 

The instability demonstrated in this letter offers a new avenue for understanding 
cosmological data, in the sense of new constraints, as well as a potential mechanism for 
generating anisotropies at specific scales and increase the tensor mode contribution and 
non-gaussianities. In this regard, much work obviously remains to be done before two- 
stream instability might be viewed as a viable, cosmological mechanism: we need to 
consider flows at arbitrary angles, include dissipation, how relative flows may develop, 
back-reaction of instabilities on the whole system, and so on. 

Finally, we note that two-stream instability is just one example of how multi-fluid 
dynamics differs from that of a single fluid. We have only considered a simple two-fluid 
model, with many features left out, but it still illustrates well the possibilities. Perhaps 
the key point is that the two-stream mechanism cannot operate in the various one-fluid 
systems that are sometimes called "multi-fluids". Although these models have several 
different constituents, they do not account for relative flows. The example considered 
here clearly demonstrates why relative flows may have interesting consequences, and 
motivates further studies of the implications. The Cosmological Principle demands a 
frame in which all constituents are at rest, but we believe strict adherence is too severe, 
may limit progress, and prevent new insights into the structure and evolution of the 
Universe. 
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Figure 1 : Time evolution of the cosmological background quantities u B and C ns (upper), and u s and C S n 
(lower), as functions of the e-fold number N for a = 1 and /} = 4/3 (matter to radiation transition). With 
these values of a and /?, the parameter k s is dimensionless while r ns has the dimensions of m 4 -3(°~n-o-s) ( a n 
functions and variables can be made dimensionless by means of a proper rescaling with some power of m). 
For all the graphs, we use *: s = 1, T ns = 0. lm 4_3 '°" n_crs ', cr n = 1.1, and <x s = 1.1. The initial values for the 
evolutions are n(0) = 3.9 x 1(T 7 , s(0) = 1, m„(0) = 0.99, h s (0) = -4.25 x 1(T 7 , A,(0) = 2, and ff,(0) = 2.89 
(for each i e [x, y, z}). 
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Figure 2: Real parts of the perturbation mode speed <x- solutions of Eq. (6), as a function of N throughout 
the transition (upper) for the same parameter values as in Fig. [T] The lower figure is a magnification of the 
region where two of the modes merge (to be distinguished from the points where lines simply cross, since 
they have zero imaginary parts, cf. figure[3j. In the limit of large N, one recovers asymptotically the four 
modes of constant speeds ± 1 / V3 and 0, as expected for radiation and matter with zero relative velocity. The 
two solutions that merge (the lower panel) pick up imaginary components (cf. figure [3) and thus correspond 
to the epoch in which the instability initiates, develops, and ends. 



8 



0.006 



0.004 



Z, 0.002 



'5b -0.002 



-0.004 



-0.006 




0.34 



0.35 



N 



Figure 3: Imaginary parts of the perturbation mode speed tr- solutions of Eq. {6), for the same parameter 
values as in Fig. [2] as a function of N throughout the epoch where the modes are unstable (ignoring the rest 
of the time evolution since all the modes are real). This corresponds to the region in Fig. |2](lower panel) 
where the real parts become degenerate: the number of parameters needed to describe the perturbations is 
therefore constant during the transition. 



9 



